Universality of the glassy transitions in the two-dimensional ± J Ising model 
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We investigate the zero-temperature glassy transitions in the square-lattice ± J Ising model, with 
bond distribution P{Jxy) = p5{Jxy — J) -I- (1 — p)S{Jxy -\- J); p = 1 and p = 1/2 correspond to the 
pure Ising model and to the Ising spin glass with symmetric bimodal distribution, respectively. We 
present finite-temperature Monte Carlo simulations at p = 4/5, which is close to the low-temperature 
paramagnetic-ferromagnetic transition line located at p ~ 0.89, and at p = 1/2. Their compari- 
son provides a strong evidence that the glassy critical behavior that occurs for 1 — po < P < Po, 
po ~ 0.897, is universal, i.e., independent of p. Moreover, we show that glassy and magnetic modes 
are not coupled at the multicritical zero-temperature point where the paramagnetic-ferromagnetic 
transition line and the T = glassy transition line meet. On the theoretical side we discuss the 
validity of finite-size scaling in glassy systems with a zero-temperature transition and a discrete 
Hamiltonian spectrum. Because of a freezing phenomenon which occurs in a finite volume at suffi- 
ciently low temperatures, the standard finite-size scaling limit in terms of TL^'" does not exist: the 
renormalization-group invariant quantity ^/L should be used instead as basic variable. 

PACS numbers: 75.10.Nr, 64.60. Fr, 75.40.Lk, 75.40.Mg 
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I. INTRODUCTION 

The ± J Ising model tl| is a standard theoretical labo- 
ratory to study the effects of quenched disorder and frus- 
tration on the critical behavior of spin systems. We con- 
sider the two-dimensional (2D) ±J Ising model defined 
on a square lattice by the Hamiltonian 



Is 
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(1) 



where ax — ±1, the sum is over all pairs of lattice nearest- 
neighbor sites, and the exchange interactions Jxy are un- 
correlated quenched random variables, taking values ± J 
with probability distribution 

PiJxy) = pSiJxy - J) + (1 - p)SiJxy + J). (2) 

For p ~ 1 we recover the standard Ising model, while 
for p = 1/2 we obtain the so-called bimodal Ising glass 
model. For p ^ 1/2, the disorder average of the cou- 
plings is given by \Jxy\ = J(^V — 1) 7^ and ferromag- 
netic (or antiferromagnetic) configurations are energeti- 
cally favored. Note that its thermodynamic behavior is 
symmetric under p — >■ \—p. In the following we set J = 1 
without loss of generality. 
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FIG. 1: Phase diagram of the square-lattice ± J Ising model 
in the T-p plane. The phase diagram is symmetric under 
p — >■ 1 — p. 



The 2D ±J Ising model has been extensively inves- 
tigated; see, e.g., Refs. 0-|j| for recent reviews. As 
sketched in Fig. [U at finite temperature it presents a 
paramagnetic and a ferromagnetic phase. They are sep- 
arated by a transition line, which starts at the pure Ising 
transition point, at p = 1 and T\s = 2/ln(l -f \/2) = 
2.26919..., and ends at the disorder-driven ferromagnetic 
T = transition, at d Q po ~ 0.897. The point where 
this transition line meets the so-called Nishimori (N) 
line [3, Si, atJlTM = 0.9527(1) and pm = 0.89083(3) 
(see also Refs. 0,13 fo'^ analytical estimates of Tm,Pm) 



is a multicritical point (MNP) fTl|,[i3- The MNP divides 
the paramagnetic-ferromagnetic (PF) transition Une in 
two parts. The PF transition line from the Ising point at 
p = 1 to the MNP is controlled by the Ising fixed point. 
Here disorder gives only rise to (universal) logarithmic 
corrections to the standard Ising critical behavior; see, 
e.g., Ref. [131 ^^^ references therein. The slightly reen- 
trant low-temperature PF transition for T < Tm belongs 
instead to a different strong-disorder (SDI) universality 
class [l,[li,[ll. 

Several studies [l6l - l28l | have discussed the possible ex- 
istence of a glassy transition, considering in most of the 
cases models with symmetric disorder distributions such 
that [Jxy] = 0. At variance with the three-dimensional 
case, no finite-temperature glassy phase occurs and a crit- 
ical behavior is only observed at T = 0. Moreover, recent 
results for the bimodal Ising model |25l - [28j have provided 
compelling evidence that the bimodal Ising model and 
other models with symmetric continuous disorder distri- 
butions, for instance the Gaussian distribution, undergo 
a zero-temperature glassy transition in the same univer- 
sality class. In all cases, for T — >■ the correlation length 
increases as T""^ with [29[ v ~ 3.55. Even though the 
asymptotic critical behavior is the same, for small values 
of the temperature and in finite volume the bimodal Ising 
glass model, and in general any model with a discrete dis- 
order distribution, presents a peculiar phenomenon re- 
lated to the discreteness of the Hamiltonian spectrum. 
This implies that, although the bimodal and Gaussian 
Ising glass have the same critical behavior in thermody- 
namic limit [2^, i.e., if one takes L -^ oo before T -> 0, 
much more care has to be exercized in the finite-size scal- 
ing (FSS) limit. In particular, as we shall see, in the bi- 
modal Ising glass model one cannot observe the standard 
FSS limit in terms of the scaling variable TL^/'^ . An ap- 
propriate variable is instead the renormalization-group 
(RG) invariant quantity (,/ L. 

In the case of the 2D ± J Ising model, a natural sce- 
nario is that a zero-temperature glassy transition occurs 
for any p in the range \—pa<p<po, and that the glassy 
critical behavior is independent of p. This implies that 
a nonzero [Jxy] is irrelevant for the critical behavior, as 
found in mean- field models [30| and in the 3D ± J Ising 
model [Hj. 

In this paper, which completes a series of papers [3, 
\\% [l3j devoted to the study of the phase diagram and 
critical behavior of the 2D ± J Ising model, we investi- 
gate the glassy behavior for 1 — po < P < Po- For this 
purpose, we present Monte Carlo (MC) simulations at 
p = 4/5, which is relatively close to the low-temperature 
PF line (po ~ Pmnp ~ 0.89), and a.t p ~ 1/2, up to 
lattice sizes L = 64 and for T > 0.1. As we shall see, 
our results provide a strong evidence of the universal- 
ity of the glassy zero-temperature critical behavior and 
thus provide strong support to the scenario of a univer- 
sal glassy critical line for T — and 1 ^ po < p < pq. 
Moreover, we provide evidence that the magnetic and 
glassy behaviors at the T = multicritical glassy point 



(MGP), where the low-temperature PF and the T — 
glassy transition lines meet, at po « 0.897, see Fig. [U 
are decoupled. Finally, we discuss the critical behavior 
of the overlap quantities along the PF line that connects 
the MNP to the T — MGP: we observe an apparently 
T-dependent critical behavior. 

The paper is organized as follows. In Sec. [ll] we de- 
fine the quantities we have considered in the MC simu- 
lations. In Sec. mil we discuss the behavior at the T = 
glassy transition. In particular, we discuss the freezing 
phenomenon observed in a finite volume at very small 
temperatures due to the discreteness of the Hamiltonian 
spectrum, the universality of the glassy critical behavior, 
and the critical behavior of the overlap susceptibility. In 
Sec. IIVI we discuss the critical behavior of the overlap 
correlations along the low-temperature paramagnetic- 
ferromagnetic transition line, below the MNP. Finally, 
In Sec.|V]we present our conclusions. 



II. DEFINITIONS 



The critical modes at the glassy transition are those 
related to the overlap variable Qx = <Tx ox , where the 

(i) 

spins (Ti belong to two independent replicas with the 
same disorder realization {Jxy}- In our MC simulations 
we measure the overlap susceptibility x and the second- 
moment correlation length ^ defined from the correlation 
function Go{x) = [{qoqx)] — [(cqctk)^], where the angu- 
lar and the square brackets indicate the thermal average 
and the quenched average over disorder, respectively. We 
define x = J2x Go{x) and 
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Go{Q)~Go{p) 



4sin^(p,„i„/2) Go[p) 



(3) 



where p = (pminjO), Pmin = 27r/L, and Go{q) is the 
Fourier transform of Go{x). We also consider some quan- 
tities that are invariant under RG transformations in the 
critical limit, which we call renormalized couplings. We 
consider the ratio S,/ L and the quartic cumulants 
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where pk = ( {Y.x1^ )^)- 

In the case of a T = transition with a nondegenerate 
ground state, as expected in the 2D Ising glass model 
with a Gaussian disorder distribution, we have X ~ ^^ 
for T = 0, hence the corresponding overlap-susceptibility 
exponent 77 vanishes, f? = 0, and 



f/4 



1, 



U2 







(5) 



for T — > 0. In particular, U22 — > indicates the self- 
averaging of the ground-state distribution, as already 
suggested by the results of Ref. 32]. Moreover, since 
77 = 0, it is natural to conjecture that the two-point over- 
lap function becomes essentially Gaussian in the limit 



T — >■ 0. If this occurs, we also have ^/L — >■ oo. As we 
shaU see, the resuhs for the ± J Ising model are consistent 
with these predictions. 

We also consider magnetic quantities. We define the 
magnetic susceptibility Xm and the second-moment cor- 
relation length ^m in terms of the magnetic two-point 
function 



Gm{x) = [(o-qctj;)]. 



(6) 



For symmetric disorder distributions we have [331 Xm = 1 
and ^,„ = for any T. For other values of p, we expect 
them to converge to a finite nonuniversal value. We also 
consider the four-point magnetic susceptibility X4„i de- 
fined by 
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(7) 
(8) 



For symmetric disorder distributions we have [33| 

X4rn = 4 - 6X. (9) 



Assuming universality we expect X4„ 
symmetric disorder distributions, i.e. 



~ X also for non- 
for any p ^ 1/2. 



III. RESULTS AT THE GLASSY TRANSITIONS 

We perform MC simulations of the square-lattice ± J 
Ising model with periodic boundary conditions for p — 
A/b and p = 1/2 and for several values of the lattice size 
L, with 8 < L < 64. We employ the Metropolis algo- 
rithm, the random-exchange method |34l |. and multispin 
coding. Furthermore, for the largest lattices we use the 
cluster algorithm described in Ref. [13|. For each lattice 
size we collect data in the range T^^-^ I^T <1.\. Atp = 
4/5 we take T,ni„ = 0.1 for L < 32, Tmin = 1/2.6 « 0.38 
for L = 48, T„ii„ = 1/2.7 « 0.37 for L = 64. At p = 1/2 
we take Tmjn = 0.1 for L < 16 and Tmjn = 1/3.3 for 
24 < L < 64. Typically, we consider 10^ disorder sam- 
ples for each T and p. In a few cases, we consider 10'^ 
disorder samples. In the following, we first discuss the 
freezing regime, which occurs for sufficiently low temper- 
atures in any finite system, then we prove the universality 
of the glassy transition by considering the FSS behavior 
of the renormalized couplings, and finally discuss the be- 
havior of the overlap susceptibility and of the magnetic 
quantities. 



A. The frozen regime 

In Figs. [H and [3] we show the MC estimates of £,/L, 
U4, and C/22- We note that the data corresponding to 
different lattice sizes cross each other around T w 0.3 and 
are mostly independent of T for T < 0.3. Similar results 
for C/4 using the bimodal distribution were also reported 
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FIG. 2: Phenomenological couplings ^/L, U4, and U22 versus 
T at p = 4/5. 



in Ref. [17]. Usually, a crossing point corresponds to a 
transition point. Instead, in the present case in which 
the disorder variables are discrete, the crossing is due 
to a nonuniversal phenomenon which is related to the 
discreteness of the Hamiltonian spectrum [25|, [2^] ■ 

To review the argument, let us consider the states cor- 
responding to the two lowest energy values for a given 
lattice size L. Their energies differ by A = i?i — i^o = 4 
and their degeneracies are given hyNo{L) and Ni{L), 
respectively. Numerical studies [21[ have shown that 
InNi/No « 41nL. At sufficiently low temperatures only 
the states with the lowest energy contribute to the ther- 
modynamics. This occurs when Nq{L) '3> Ni{L)e^^/'^ , 
i.e., for 
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FIG. 3: Phenomenological couplings £,/L, U4, and U22 versus 
T at p = 1/2. 



In this regime the observed behavior is independent of T. 
In the opposite Umit, i.e., when Nq{L) <^ Ni{L)e~'^''^ , 
the presence of the gap is negligible and the system is 
expected to have the same behavior as models with con- 
tinuous distributions. The crossover from one regime to 
the other occurs at a L-dependent freezing temperature 
Tf{L) which scales as 1/ InL. It is natural to define Tf{L) 
by requiring No{L) = Ni{L)e~^/'^f , but this definition 
is somewhat unpractical. In practice, Tf(L) can be esti- 
mated from the data by identifying it with the tempera- 
ture that marks the onset of the T-independent behavior 
of the different observables. For p — 4/5, the estimates 
of ^/L and C/4 reported in Fig. [2] allow us to estimate 
Tf{L) « 0.41, 0.36, 0.32, 0.28 for L = 8, 12, 16, 24, re- 
spectively. Slightly larger results are obtained by using 
1/22 ■ The estimates of Tf{L) for p = 1/2 are close to 
those obtained for p = 4/5, showing that Tf{L) is little 



dependent on p. Consistently with the above- reported ar- 
gument, the freezing temperature Tf[L) approximately 
decreases as 1/lnL, see Fig. SI A fit to c/lnL gives 
c ss 0.9. In the frozen region the estimates of the renor- 
malized couplings should be very close to the correspond- 
ing T = estimates, since they are essentially determined 
by the lowest-energy configurations. The data are con- 
sistent with this prediction. Indeed, see Fig. |4l for both 
p = 4/5 and p = 1/2, Ua and U22 below T/ slowly ap- 
proach the values U22 — and C/4 = 1, respectively, 
for L — > 00. In particular, U22 apparently vanishes as 
U22 ^ 1/ Ini. Below T/ the ratio S,/L increases as L, see 
Fig. [5l indicating that ^ ~ L^ at T = 0. These results 
imply that, at T = 0, the large-L limit of the renormal- 
ized couplings is identical to that observed in models with 
continuous distributions; in this case, as we discussed be- 
fore, we predict U22 ^- 0, C/4 ^- 1, and ^/L — >■ 00. This 
equality should not be taken as an obvious fact. For in- 
stance, the stiffness exponent is different in the two cases. 



B. Finite-size scaling in the presence of freezing 

The presence of freezing for T < Tf{L) makes the 
study of the T = glassy critical behavior quite hard. 
Indeed, in order to observe the glassy critical behavior in 
Ising glass models with a discrete distribution, one must 
approach T = by keeping T ^ Tf{L) for each lattice 
size. This makes a standard FSS analysis impossible. In- 
deed, in the FSS limit a RG invariant quantity R should 
scale as 



R = fniTL 



1/l.N 



(11) 



The condition T ^ Tf{L) implies that this scaling be- 
havior can only observed for 



y^i/z. ^ Tfi^L^'" 



InL ' 



(12) 



For L — )■ oo, the ratio L^/'^/lnL diverges and thus this 
makes the range of values of TL^^'^ which are accessible 
smaller and smaller as L increases. This implies that the 
standard FSS limit, T ^ 0, L ^ oo at fixed TL^/" does 
not exist. However, as we shall now discuss, one can still 
study FSS if one uses the ratio £^/L as basic FSS variable, 
i.e., if one considers the scaling form 



R = gR{^/L). 



(13) 



Usually, expressions ([TTt and (IT^ are equivalent. This is 
not the case here: only the FSS scaling form ([13]) holds 
in the presence of freezing. As is clear from Figs. [2] and [3l 
the ratio ^/L at fixed L increases as T decreases. Hence, 
the condition T ^ Tf{L) translates into 






(14) 



where ^/ is the value of ^ in the frozen region. For L ^ oo 
^f/L diverges and thus, by increasing L, one has access 
to the whole FSS region. Thus, in the presence of freezing 
FSS can still be used but only in the form (IT^ . Since 
v does not appear in Eq. ^3]) . FSS cannot be used to 
determine v (this exponent can only be determined by 
using infinite- volume data), though it can still be used 
to check universality. Note that the presence of freezing 
and the limitations in the use of FSS are always expected 
in models with a T — transition and discrete Hamil- 
tonian spectrum. In particular, these phenomena should 
also be considered in the 3D diluted ± J Ising model close 
to the percolation point, where the glassy transition tem- 
perature vanishes [35i] . 



C. Renormalized couplings and universality 
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In order to verify universality, we consider the quartic 
cumulants t/4 and f/22j which should scale according to 
Eq. (J13p . The function gn^x) should be universal, hence 
p-independent. Universality is nicely supported by the 



FIG. 7: Renormalized U22 versus ^/L for p = 4/5. We plot 
data for all available values of L. Below, we only show the 
results for ^/L > 0.3. 



data shown in Fig. |6l As expected, we find that U4 —> 3 
and U22 -> for £,/L -> 0, and t/4 -J> 1 and U22 -^ for 
^/L -^ cx). All data for C/4 fall onto a single curve with 
small scaling corrections, supporting the existence of the 
finite-size scaling limit in terms of £,/L, as discussed in 
the previous section. The convergence to a single curve 
is also clear in the case of U22, although corrections are 
quite evident. Note that, in the region around the peak, 
for ^/L « 0.3, the data at p = 4/5 and p = 1/2 converge 
from opposite sides. In these universality checks there are 
no free parameters to be adjusted, and thus these com- 
parisons provide strong support to the hypothesis that 
these models belong to the same universality class. An 
analogous universality check was performed in Ref. [2^. 
By studying the FSS behavior of U^ versus (,/L, it pro- 
vided numerical evidence that the ± J Ising model with a 
bimodal distribution (our model with p — 1/2) and that 
with a Gaussian distribution of the couplings belong to 
the same universality class. It is interesting to note that 
while the data for C/4 scale nicely up to ^/L < 0.9, for U22 
significant deviations are observed at smaller values of 
^/L, see Fig. [71 For p — 4/5 and L = 8 the data show sig- 
nificant deviations close to the peak, then approach the 
common curve and then show again a significant devia- 
tion — the data turn up — for ^/i > (^/i)max ~ 0.55. A 
similar phenomenon occurs for L — 12. For L — 16 devia- 
tions close to the peak arc quite small, but again the data 
begin to turn up as £^/L > (^/i)max ~ 0.6. For L = 24 
FSS holds quite nicely, at least up to (^/-/j)max ~ 0.65. 
The value (^/i)max marks the onset of the crossover re- 
gion between the critical regime where FSS holds and 
the freezing regime that sets in at £,f/L. Note that £,f/L 
is significantly larger than {S,/L)^ax, indicating that the 
breaking of FSS occurs before freezing. For p — 1/2 
the conclusions are similar, although, for a given L, 
(^/Z/)inax is significantly smaller than the corresponding 
value for p = 4/5. For instance, for L = 16, we have 
(e/i)max « 0.45 for p = 1/2 and (C/i.)max ~ 0.60 for 
p = 4/5. This clearly reflects the fact that ^/ for p = 1/2 
is smaller than for p — 4/5, see Fig. [5j 



D. Overlap susceptibility and exponent -q 

Finally, we investigated the critical behavior of the 
overlap susceptibility. As discussed in Ref. [31|] it should 
behave in the FSS limit as 



X^uliT)L'-^F^{aL), 



(15) 



where Uh is an analytic function of T which is related to 
the overlap-magnetic scaling held. In order to determine 
r] we have performed fits of the data with p = 4/5 to 



Inx = (2 - 77) Ini + PniT) + QmiC/L), 



(16) 



where P„(a;) and Qmi^) Siie polynomials in x with 
-Fn(O) = 0. To avoid any bias from the presence of the 
freezing region, we have only used the data satisfying 



C/L < (^/L)max, where, for L < 24, (C/i)max is the value 
determined before from the analysis of 1/22- For L — 32 
we used somewhat arbitrarily (^/£)max = 0.7, while for 
L = 48, 64 we used all our data which in any case sat- 
isfy ^/i < 0.65. To identify scaling corrections we only 
considered data satisfying T < Tmax and L > Lmin for 
several values of Tmax and Lmin- The results depend 
strongly on these parameters. For Tmax = 1-2 we obtain 
7] = 0.39(1), 0.33(3) for L^in = 8, 16. For £„,;„ = 16, we 
obtain 77 = 0.27(3), 0.24(3), 0.22(4) for T^^y, = 1,0.8,0.6. 
Apparently the results always decrease as Tmax decreases 
and imin increases. It is impossible to estimate reliably 
77 from these results and thus we only quote an upper 
bound: 



^^0.2 



(17) 



In principle, the same analysis can be performed for p = 
1/2. However, in this case the estimates of {^/L)^^^ are 
smaller, so that fewer data can be used in the fit. In 
practice, no estimates of 77 can be obtained. 

As we mentioned at the beginning there are strong the- 
oretical reasons to expect 77 = 0. Thus, we tried to verify 
whether our results are consistent with this hypothesis. 
For this purpose we considered the data with Tmax = 1.0 
and imin = 12 and we fitted them to Eq. P^ setting 
7; = 0. We find that the data satisfying ^/L < (<^/i)max 
are reasonably described by Ansatz P^ . In Fig. [5] we 
report 



Xr 



XL- 



^PAT) 



(18) 



where Pn(T) is the polynomial determined in the fit 
P^ . The agreement is quite good up to {^/L) w 0.65. 
It should be noted that, if we also include data with 
^/L > (C/i)max, the fits become much less dependent 
on Tmax and imin, excludc 77 = 0, and give the esti- 
mate 77 w 0.2: Apparently the data that belong to the 
region (i^/i)max < £./L < £.f/L are well described by 
Eq. P^ with 77 = 0.2, while they cannot be fitted by 
taking 77 = 0. This is evident from Fig. [5] where we 
report x-^~^'^exp[— P„(r)]: in this figure essentially all 
data fall on a single rescaled curve. Note that the quality 
of the collapse is only marginally better than that given 
in Fig. [5]— the xVDOF of the fit (DOF is the number 
of degrees of freedom of the fit) is similar in the two cases 
— although this is not evident from the figures since they 
have a completely different vertical scale. 

As a final check we verify if the data at p = 0.5 are 
also consistent with 77 = 0. The results are reported in 
Fig. [TUl Also in this case the data with ^/L < 0.65 are 
consistent with 77 = 0. In the plot we have also multiplied 
Xrosc by a constant in such a way that Xrcsc assumes the 
value Xresc ~ 56 for ^/L « 0.3 as it does for p — 4/5. 
With this choice the curves for Xrosc should be the same 
for both values of p. As it can be seen the shape of the 
two curves is indeed the same. Quantitatively, the two 
curves are the same up to £,/L sa 0.45, while they differ 
significantly for ^/L w 0.6, 0.7. This is not surprising. 




FIG. 8: Rescaled susceptibility Xrcsc defined in Eq. (|18[) for 
p = 4/5 and i] — 0. 
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FIG. 11: The ratio g™/x at p = 4/5. 




FIG. 9: Rescaled susceptibility Xrcsc defined in Eq. (|18p for 
p = 4/5 and r? = 0.2. 



As we have explained, for p = 1/2, the data such that 
^/L ^0.5 are probably already in the crossover region 
before the onset of freezing. 



E. Magnetic quantities 

The magnetic variables do not become critical in limit 
T — > 0. Indeed, the magnetic susceptibility Xm and sec- 
ond moment correlation length ^^ are finite in the limit 




FIG. 10: Rescaled susceptibility Xrcsc defined in Eq. (|18|) for 
p = 1/2 and r? = 0. 



T ^- 0. For p ~ 1/2 we have Xm = 1 and ^^ = for any 
T. For other values of p, they converge to nonuniversal 
values such that Xm > 1 and £,m > 0. For p = 4/5 we 
find Xm ~ 24 and 6™ « 3.0 in the limit T ^ 0. We 
also consider the four-point magnetic susceptibility X4m 
which should scale as x in the critical limit. In Fig.fTTJwe 
plot the ratio g^/x, where gm = -X4m/(XmCn)- This 
quantity shows smaller scaling corrections than Xim/x 
and clearly converges to an L independent constant. The 
asymptotic behavior sets in for T < 0.5, before the freez- 
ing region T < 0.35. 



IV. THE OVERLAP CRITICAL BEHAVIOR 

ALONG THE LOW-TEMPERATURE 

PARAMAGNETIC-FERROMAGNETIC 

TRANSITION LINE 

We now investigate the behavior of the overlap cor- 
relations along the low-temperature PF transition line, 
see Fig. [H from the MNP to the T = axis. The crit- 
ical behavior of the magnetic correlations was numeri- 
cally studied in Ref. [4] by varying p for two values of T: 
T = 1/1.55 and T = 1/2. It was found that the criti- 
cal behavior was universal, controlled by a single strong- 
disorder fixed point, with critical exponents v sa 3/2 and 

Tim ~ 1/8. 

The ferromagnetic T = transition point at p = po ~ 
0.897, where the low-temperature PF transition line ends, 
is a multicritical glassy point (MGP), because it is con- 
nected to three phases and it is the intersection of two 
different transition lines, the PF line at T > and the 
glassy line at T = 0. At T = the critical point at 
p = Po separates a ferromagnetic phase from a T = 
glassy phase, while for T > the transition line separates 
a ferromagnetic from a paramagnetic phase. Therefore, 
on general grounds, the critical behavior at T = and 
p = Po should differ both from that observed along the 
PF line and that observed along the glassy line T = 0, 
p > Po, unless the magnetic and glassy critical modes are 
effectively decoupled. Such a decoupling is appare ntly 
supported by the numerical results of Refs. (a. la. uA. l36| . 



Indeed, the estimates of magnetic critical exponents at 
T = are quite close and substantially consistent with 
those found along the PF transition line below the MNP. 
All results are therefore consistent with a single magnetic 
fixed point that controls the magnetic critical behavior 
both at T > and at T = 0. The analysis of the over- 
lap correlation functions also supports the decoupling of 
the critical modes. Indeed, the arguments we gave in 
Sec. In] on the behavior of the (overlap) renormalized cou- 
plings and correlation functions should hold at T = for 
any value of p — hence, in the ferromagnetic phase, at 
the MGP as well as along the glassy transition line — 
since they only rely on the assumption of a nondegen- 
erate ground state. Therefore, also at the T = MGP 
we expect C/4 = 1, U22 = 0, ^/L = 00, and 77 = 0. The 
behavior of the overlap correlations is therefore identical 
for p > po, for p < pq, and at the MGP, in agreement 
with the decoupling scenario. 

We wish now to understand the critical behavior of 
the overlap correlations along the PF line. In order 
to investigate this issue, we perform MC simulations at 
two critical points along the low-temperature PF line, at 
[T = 1/1.55 « 0.645, p = 0.8915(2)] and [T = 1/2, p = 
0.8925(1)], as determined in Ref. ^], for lattice sizes up to 
L — 48. The FSS analysis of the data of the overlap sus- 
ceptibility X shows that x ~ L^~^ with quite small but 
nonzero values of rj. Fits of the data satisfying L > 16 
which take into account the nonanalytic scaling correc- 
tions give ri = 0.046(6) at T = 1/1.55 and -q = 0.038(4) 
at T = 1/2. Note that these values are much smaller 
than the pure Ising value rj = 1/2 (which is simply twice 
the value of the magnetic exponent rjm = 1/4) holding 
along the PF line from the pure Ising point to the MNP, 
and also much smaller than the value 77 = rjm — 0.177(2) 
at the MNP point. The FSS analyses of the renormal- 
ized couplings ^/L, U4 and U22 lead to apparently T- 
dependent critical values: U22 = 0.018(2), U4 = 1.044(4), 
£./L = 1.7(1) at T = 1/1.55, and U22 = 0.008(2), 
C/4 = 1.025(2), £,/L = 2.2(1) at T = 1/2 (where the 
errors are essentially due to the uncertainty on pc and on 
the scaling correction exponent). Note that these values 
are very close to the values at the T = MGP; still we 
consider unlikely that the overlap behavior is the same as 
that occurring at the MGP, mainly because this would re- 
quire S^/L = 00 along the whole PF line between the MNP 
and the MGP. Indeed, the condition ^/L = 00 is quite un- 
likely for a finite-temperature transition. These results 
can be better explained by a T-dependent asymptotic 
critical behavior which, with decreasing T, approaches 
the T — glassy behavior characterized by the values 



77 = 0, U22 = 0, t/4 - 1, ^L = C30. 



V. CONCLUSIONS 

In this paper we consider the two-dimensional ± J Ising 
model, focusing mainly on the T — glassy transition 
occurring for 1 — po < P < Po- The main results are the 
following. 

i) We first discuss the freezing phenomenon that oc- 
curs on any finite lattice at sufficiently low tem- 
peratures. We investigate the behavior of several 
quantities in this regime, verifying explicitly the 
expected logarithmic dependence on the lattice size 
L. We also show that the presence of this regime 
makes it impossible to use the standard form of 
FSS: the FSS fimit T ^ 0, L ^ 00 at fixed TL^^" 
does not exist. FSS can be formulated only if one 
considers £,/L as basic FSS variable. 

ii) We study the FSS behavior of the renormalized 
couplings U4 and U22 for p = 4/5 and p = 1/2 
as a function of ^/L. We find that they have the 
same FSS curves, a clear indication that the critical 
behavior for p = 4/5 and p = 1/2 is the same. This 
allows us to conjecture that the critical behavior is 
independent of p in the interval 1 — po < p < Po , 
and therefore, together with the results of Ref. j26| . 
that there exists a single universality class for 2D 
Ising glassy transitions. We also investigate in de- 
tail the critical behavior of the overlap susceptibil- 
ity, showing that the numerical data are consistent 
with ?7 = 0, if one discards data that are close to 
the region where freezing occurs. 

iii) Finally, we discuss the critical behavior of the over- 
lap variables along the PF line. An analysis of the 
numerical data available for the T — MGP indi- 
cates that at this point glassy and magnetic modes 
are decoupled. For T > we observe an apparent 
T-dependent critical behavior. Note that a similar 
phenomenon was also observed in the XY model 
with random shifts 37]: along the critical line that 
starts at the XY pure point and ends at the Nishi- 
mori multicritical point, magnetic quantities show 
a universal behavior while overlap variables show a 
disorder-dependent critical behavior. 
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